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Abstract Isotherm curves of water sorption for wood at 
various temperatures were analyzed based on the dual 
mode theory where the total coverage was represented by a 
linear combination of the Langmuir and Henry equations. 
The saturation concentration and affinity constant of the 
Langmuir equation and the parameter of Henry’s law had a 
transition point near 60°C. The analysis based on the dual 
mode theory found that the constants for whole wood were 
related to those of wood components and depended 
more on their glass transition temperatures. That is, it was 
theoretically demonstrated that the characteristic tem¬ 
perature dependence of water sorption for wood occurs 
because wood consists of three components (cellulose, 
hemicellulose, and lignin) with different glass transition 
temperatures. 
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Introduction 

The isotherm curve of water sorption for wood is sigmoid 
and similar to other natural products. 13 Skaar 4 introduced 
various sorption theories and described how the water sorp¬ 
tion of wood fit the Hailwood and Horrobin theory, 5 which 
is a kind of dual mode theory. This dual mode theory has 
been applied to both natural and chemically modified 
wood. “ The general dual mode theory, which was pro¬ 
posed by Meares 11,12 and Michaels et al., 13 has been applied 
successfully to the gas adsorption of synthetic polymers, 1417 
in which the total coverage is represented by a linear com¬ 
bination of the Langmuir and Flory-Huggins equations. In 
some cases, Henry’s equation is exchanged for the Flory- 
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Huggins equation because the former is the approximate 
equation of the latter. 18 This theory differs from the 
Hailwood and Horrobin theory in that the Langmuir 
sorption process is dependent on the Henry sorption pro¬ 
cess for the latter theory. 

Nevertheless, the characteristic temperature depen¬ 
dence of water sorption for wood has not been examined in 
detail yet based on the dual mode theory, although it has 
been applied to various materials. Therefore, this report 
analyzes the temperature dependence of water sorption 
for wood based on the dual mode theory considering 
wood components that have different glass transition 
temperatures. 


Data used in the analysis 

The isotherm data for 0°-100°C used in this analysis were 
taken from the Wood Handbook published by the US 

i o 

Department of Agriculture and are probably the standard 
isotherm data for wood. Figure 1 shows typical isotherm 
curves from the data. 


Results and discussion 

Application of the dual mode theory and temperature 
dependence of its parameters 

In the dual mode theory, total moisture content (cm 3 STP/ 

A "1 "1 1 O 

cm’) is represented by the following equation: “ J 

c = c D +c H (1) 

where C D and C H are described by the Henry and Langmuir 
equations, respectively, that is 

C D = K D h (2) 

Chb'h 
H 1 + b'h 


( 3 ) 
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Fig. 1. Typical sorption isotherm curves at various temperatures using 
data from the Wood Handbook published by the US Department of 
Agriculture. 19 m, Moisture content, h, relative humidity 



Fig. 2. Relationships between the ratios of relative humidity to mois¬ 
ture content {him) and relative humidity ( h) 


where h is the fractional relative humidity, K D is Henry’s 
parameter, C' H and b 'are the saturation concentration and 
modified affinity constant in the Langmuir equation, respec¬ 
tively, such that b'— bp 0 , where b and p 0 are affinity constant 
and vapor pressure, respectively. When using the moisture 
content units (g/g) generally used for wood, 

m = m D + m H (4) 

with 



K D h 
1 -K D h 


m H = 


mjo'h 
1 + b'h 


( 5 ) 

( 6 ) 


where k is a constant by unit conversion, and ra 0 is the 
saturation concentration (g/g). Equations 5 and 6 show the 
linearity of h/m D vs h and h/m H vs h. Note that the dual 
mode theory differs from the Hailwood and Horrobin 
theory in that ra H and m D are dealt with independently. In 
the Hailwood and Horrobin theory, the constants ra 0 , b\ 
and K d are calculated by solving a quadratic equation ob¬ 
tained by substituting Eqs. 5 and 6 for Eq. 4. 

Figure 2 shows him vs h for wood from the Wood Hand¬ 
book published by the US Department of Agriculture 19 and 
that the plots are not symmetric. This implies that the 
Hailwood and Horrobin theory cannot be applied to the 
curves in Fig. 2, because him vs h must be a quadratic curve 
in the Hailwood and Horrobin theory. Nevertheless, in Fig. 
2, these curves are linear in both the lower and higher 
relative humidity regions, i.e., they are approximated by the 
Langmuir theory at low relative humidity and by Henry’s 
law at high relative humidity. As a result, we can assume 
that isotherms at various temperatures for wood are de¬ 
scribed by a linear combination of the Langmuir theory and 
Henry’s law. The isotherm obtained using Eqs. 4-6 is in 
good agreement with the experimental results without the 



Fig. 3. Comparison of the experimental results (open circles) and theo¬ 
retical calculation (solid line) for a sorption isotherm curve 


intermediate region of relative humidity, as shown in Fig. 3, 
which is example one of the all data. Note, however, that 
this result is not necessarily valid theoretically, because 
wood is a multicomponent system. We must consider the 
effects of the different components. 

The temperature dependence of parameters ra 0 , b ' and 
K d is shown in Figs. 4 to 6, which were calculated from h/m H 
vs h for low relative humidity and from h/m D vs h for high 
relative humidity. We note that the slope of ra 0 vs T is 
positive because the slope is generally negative for many 
polymers. Moreover, they have a transition point near 60°C 
(= 0.003 K 1 ). This characteristic temperature dependence 
suggests that the overall water sorption of wood is affected 
by the glass transition temperature of its components, 
considering the glass transition effect on gas sorption for 
various polymers. 1517 Therefore, we will try to explain the 
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temperature dependence, that is, the appearance of the 
temperature transition seen in Figs. 4-6, regarding a multi- 
component system whose components have different glass 
transition temperatures. 

Dual mode theory and glass transition temperature 

It is well known that gas sorption behavior changes abruptly 
at the glass transition temperature T g for synthetic poly¬ 
mers. For glassy polymers below T g , both gas sorption 
on the surface in the unrelaxed volume, i.e., the sum of 
microvoids, and dissolution into the matrix proceed during 
sorption; that is, a sorption isotherm consists of contribu¬ 
tions from both Langmuir sorption and Henry dissolution. 
The unrelaxed volume below T g decreases with increasing 
temperature, so that the coverage C' H , that corresponds to 
ra 0 , decreases with temperature according to the following 
equation: 17 

C^=yAa{T g -T) (7) 

where yis a constant, and Aai s the difference between the 
cubic thermal expansion constants above and below T. 
Moreover, T g also decreases with increasing sorbate con¬ 
centration, because the flexibility of molecular chains in the 
matrix increases with the sorbate concentration. As a result, 
Eq.7 is reduced to the following equation: 

Q = YAa(T ia - PC*-T) = yAcc[T g (C*) - T ] (8) 

where T go is T g at C* = 0, and (3 and C* are a constant and the 
effective sorbate concentration, respectively. Equation 8 is 
valid for many synthetic polymers. 17 This equation shows 
that the Langmuir sorption at T < T g is affected by both the 
measurement temperature and the sorbate concentration. 
The result shown in Fig. 4 is in conflict with Eqs. 7 and 8 
because the slope of C' H vs T of equations is negative. This 
conflict is discussed hereinafter. 



-10 30 70 110 

Temperature (°C) 

Fig. 4. Temperature dependence of the saturation moisture content, 

m () 


For simplicity, we use Eq. 7, because the effect of C* on 
T g is slight at low relative humidity. By contrast, for rubbery 
polymers above T, most microvoids disappear, so that the 
sorbate dissolves in the matrix, and its dissolution obeys 
Henry’s equation or the Flory-Huggins equation in some 
cases. That is, the sorption behaviors above and below the 
glass transition temperature differ because the flexibility of 
the molecular chains changes. 

The glass transition temperatures of wood components 
differ. Many researchers have reported the glass transition 
of wood and its components under various wet condi¬ 
tions. 2027 The reported values differ because of differences 
in the measurement methods. Furuta et al. 27 reported -40°C 
for hemicellulose, 50-100°C for lignin, and above 100°C for 
cellulose using viscoelastic measurements under wet condi¬ 
tions. These values are much lower than those of Goring 20 
determined with the softening point apparatus or the me- 
chanical measurements of Salmen and Back. ’ Hillis and 

r\ a 

Rozsa" found two peaks below 100°C in the temperature 
dependence of the rotation rate under a constant torque for 
wet wood, while Irvine 25 and Kelley et al. 26 reported glass 
transition temperatures of 60°-70°C for milled wood lignin 
(MWL). Considering these results, the glass transition 
temperatures of lignin and hemicellulose are about 60°C 
and below room temperature, respectively. Therefore, the 
temperature dependence of ra 0 , b ' and K D for whole wood 
is expected to be apparent. 

Whole wood as a multicomponent system can be ana¬ 
lyzed using dual mode theory to clarify the characteristic 
temperature dependence of ra 0 , b\ and K D , as shown in 
Figs. 4-6. We will use the units cm 3 STP/cm 3 and Eq. 8 in the 
following discussion. The temperature dependence of the 
saturation concentration, affinity, and Henry’s parameter 
does not change on unit conversion. 

Temperature dependence of Langmuir’s parameters 

First, we examine isotherms for low relative humidity, for 
which the isotherms should obey the Langmuir equation, as 
shown in Fig. 2. 



1/7 (K" 1 ) 

Fig. 5. Relationship between the logarithmic values of the affinity 
constant h'and the reciprocal of absolute temperature 1 IT 
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Fig. 6. Temperature dependence of Henry’s law parameter K B 


In our discussion, the wood components hemicellulose, 
lignin, and cellulose are represented by the subscripts 1, 2, 
and 3, respectively. For example, their respective glass tran¬ 
sition temperatures are T gi , T %1 , and T g3 , and T g{ < T gi < T gy 
Assuming that each component obeys Eq. 1, then Eq. 1 for 
the entire system of wood is represented by 

C — =I>(C Df + C H (9) 

i =1 i =1 

where 0, is the volume fraction of component /. 

According to dual mode theory, sorption below T g 
causes both Langmuir sorption and Henry dissolution. 
However, the contribution of the Henry dissolution is small 
at low relative humidity, so that the sorption is approxi¬ 
mated by the Langmuir sorption. By contrast, the 
Langmuir-type sorption levels off at high relative humidity, 
so that its contribution is constant, C' H . Then, the moisture 
content of component i is represented by 

c,=(^ D ,+c;,)/i (rsr g ) (io) 

where K D . is Henry’s parameter for component i. That is, 
the sorption for high relative humidity is approximated by 
Henry’s equation. Of course, the moisture content at T ^ T 0 
obeys Henry’s law for all values of relative humidity. 

For low relative humidity, the total moisture content is 
represented by Eqs. 11-13 in each temperature region using 
the previous approximation. 



C = (j) l K Di h + '^(pi 


1=2 


Ckb'h 
l + b'h 





C = Y,<j>iK Di h + (l) 

i=2 

(l + b'M) 


C' H b’Ji 


1 + b 3 h 
(4 T 






where C' H , b' n and K D . are the saturation moisture content, 
modified affinity constant, and Henry’s parameter, respec¬ 
tively. Equations 11,12, and 13 are derived using h 1 and 
assuming b\h 2 1. In general, b\< 1 for many synthetic 
polymers. 

Note that from Eqs. 11-13, the Langmuir equation is not 
exactly valid, but is approximately valid in the region of low 
relative humidity for wood when the T g s of the components 
differ. That is, the Langmuir-type isotherm for low relative 
humidity in Fig. 2 is apparent. The constants C' H (or ra 0 ), b\ 
and K d obtained from the isotherm for whole wood reflect 
those of the components, so that their temperature depen¬ 
dence depends on the glass transition temperatures of the 
components. 

This is why the Langmuir curve in the dual mode theory 
does not level off for wood, 28 while the curve does level off 

z' 1—1 

in the Hailwood and Horrobin theory. ’ The former should 
be valid because one of the components, moist hemicellu¬ 
lose, does not obey the Langmuir equation, but obeys 
Henry’s law at room temperature (T g ^ T ^ T g2 ), as shown 
in Eq. 12, whose contribution is represented by 0 X K D{ h in 
Eq. 12. 

Next, we examine the temperature dependence of C' H 
and b' for the apparent Langmuir-type isotherm, which 
is found in the region of low relative humidity. From Eqs. 
11-13, the apparent C' H is represented by 




(15) 


C' - 

'-'H — 


exp[(zi// H . 

i =1 

(/> 3 y 3 Aa 3 b 3 T g3 jb 3 


AH bi )/RT]-(/) 3 y 3 Acc 3 T 

(4 = T = Tg 3 ) 




b{ 


J 


(16) 


where AH h . and AH D . are the enthalpy of the Langmuir 
sorption and the Henry dissolution for component i, respec¬ 
tively, and AH h . AH d . < 0. In the derivation, C' n . = yAa^T,, 
- T), b t = b ( !txpi-AH H JRTl and K D . = K^exp [~AH D JRT\ are 
used, where is a constant, R is the gas constant, and T is 
the absolute temperature. 

Equation 14 as a function of T represents a line with a 
negative slope at T < T g . Because we are concerned with 
the appearance of the transition at 60°C, as shown in 
Fig. 4, we do not deal with Eq. 14 here because T g room 
temperature. 
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The slopes of Eqs. 15 and 16 depend on the contribution 
of the exponential term in the equations, that is, the contri¬ 
bution of Henry’s dissolution. Given that exp [(AH h . - 
AH Di )/RT\ is a monotonically increasing function of T be¬ 
cause AH h . AH d . < 0, Eqs. 15 and 16 are monotonic 
functions of T and are approximated by straight lines in the 
narrow temperature region. Their slopes above and below 
r g2 differ. If the exponential term in Eqs. 15 and 16 is 
greater than the term of T, the slopes above and below T g2 
(= 60°C) are positive, and the former is less than the latter. 
This explains the relationship between ra 0 and T in Fig. 4, 
where ra 0 corresponds to C ' H . 

The above discussion explains partly the conflict be¬ 
tween Eqs. 7 and 8 and the result shown in Fig. 4. Unfortu¬ 
nately, we cannot clarify that all because the values of most 
of the constants in Eqs. 15 and 16 is not known. Neverthe¬ 
less, our discussion explains the appearance of the transi¬ 
tion at 60°C in Fig. 4 at least. 

Next, we discuss the temperature dependence of b \ The 
apparent b'tox the whole system is obtained from Eqs. 11- 
13. Considering b' = exp [~AH U JRT] in the Langmuir 
theory, b' in an each temperature region is represented by 

b' = ±b'=±b°^ V [-AH n jRT] (T § 4 ) (17) 

i =1 i =1 

b ' = t b ' = X b < exp[-^H, /RT] (4 sr§4) ( 18 ) 

i=2 i =2 

b ' = b ’ = b 3 exp[-4 // H3 /RT] (4 STSrJ (19) 

The affinity b' as a function of IIT is a monotonically in¬ 
creasing function considering AH h . < 0 so that \n(b / ) is also 
monotonically increasing. Then, we derive the following 
equations as a slope when \n(b') is differentiated by 1 IT. 


d\n(b') 

d(l/T) 











( 22 ) 

where AH Rn and AH H3 are the enthalpy of water sorption for 
lignin and cellulose, respectively. The affinity of water for 
cellulose is greater than that for lignin, and water sorption is 
an exothermic process, AH H3 < AH H2 < 0. Thus, comparing 
Eqs. 21 and 22, we have 


0 < 


dln(b') 

d{VT) 


(l/T Sl )<{l/T) 


dln(b') 

d{VT) 


(1/T)z(l/T g2 ) 



Equation 23 reveals that the slope at 1/T g2 < 1 IT is less than 
that at 1 IT < 1/T g2 for \n(b') vs 1 IT. This equation explains 
the transition at 1 IT = 0.003 K 1 (= 60°C) in Fig. 6 , if T g2 = 
60°C, which is the glass transition temperature for MWL 

s' 

according to Irvine and Kelley et al. 

Therefore, we conclude that Figs. 4 and 5 reflect 
the temperature dependence of C' H (or ra 0 ) and b' for 
wood components whose glass transition temperatures 
differ. 


Temperature dependence of Henry’s parameter 


Even in the region of high relative humidity, the isotherm 
obeys the Langmuir sorption at T < T g and the Henry 
dissolution at T g < T. However, the Langmuir sorption 
apparently levels off at high relative humidity. Therefore, 
the moisture content at T < T g is represented by Eq. 10, 
as mentioned above. Of course, it is represented by Henry’s 
equation at T g < T. Accordingly, the moisture content of the 
entire system is represented by 



Equations 24-26 show that Henry’s law apparently holds 
for the region of high relative humidity and that Henry’s 
parameter K D for the entire system reflects the sorption and 
dissolution properties of wood components. As a result, K D 
is represented by 


* D = 2>(*d, +q.)=x^ d ,+X c h, (^4) (2?) 


i =1 


i =1 


i =1 


k d =< p t K Vh + 5>,(* Dl + q,)=5> ( * d ,+ yc,', 

i=2 i =1 i=2 

(4 S T s 4) 


(28) 
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Here, the temperature dependence of C' H and K D . are 
represented by C' H = yAa(T g - 7) presented in Eq. 7 and 
= T^.exp [~AH d JRT ], respectively. 

Now, we examine K D as a function of T above and below 
r g . From Eqs. 27-29, the slope dK^/dT of 7£ D vs T is derived 
as follows: 


values of most of the constants in the derived equation are 
not known for wood. Therefore, a strict discussion of this is 
expected after these values are clarified. 


References 


= i(0,<4W D; / RT 2 )exp[-AH 0i / RT] - 

d l i =i i =i 

(30) 

/i?r 2 )exp[-zl// Di /i?T] - 

U1 i =1 i=2 

(4 = T Sk 4 ) 

( 31 ) 

^ /i?4)exp[-zi// D . /i?T] - (piYjAcijbj 

a l i=] 
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As iC D . = 0 / i^D.(zi// D ./i?r 2 )exp[-Zl// D ./i?r] is a monotonically 
decreasing function of 7, from Eqs. 31 and 32 we find 


dK D 

dT 
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< 
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<T 


<0 



Equation 33 shows the appearance of the transition at T g2 
for K d vs T. As K d . = ^K° Di (AH D JRT 2 )cxp[-AH D JRT] is 
approximated by a straight line with a negative slope in the 
narrow region of 7, K D vs T in 7 g ^ 7 ^ 7„ 3 consists of two 
straight lines with a transition point at 7 g2< This is in agree¬ 
ment with the result shown in Fig. 6. Therefore, the K D for 
whole wood depends on the thermal properties of its 
components. 


Conclusions 

The isotherm curves of water sorption at various tempera¬ 
tures for wood were analyzed using the dual mode model. 
The saturation concentration and affinity constant of 
Langmuir’s equation and Henry’s parameter as a function 
of temperature had a transition point at 60°C, which is the 
glass transition temperature of lignin. The analysis consid¬ 
ering a multicomponent system using the dual mode model 
found that the constants for whole wood are related to 
those of its components and depend on their glass transition 
temperatures. That is, the characteristic temperature de¬ 
pendence of water sorption for wood occurs because wood 
consists of components with different glass transition 
temperatures. 

In this report, the slopes of C' H and K D near and below 
7„ 2 could not be discussed in detail because the change in 
the microvoid volume plasticity with water sorption and the 
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